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FOREWORD 


The work described here was done at the Atomics International 
Division of Rockwell International Corporation, under the direction 
of the Space Nuclear Systems Division, a joint AEC -NASA office. 
Project management was provided by NASA-Lewis Research Center 
and the AEC-SNAP Project Office. 


DISTRIBUTION 


This report has been distributed according to the category 
"Systems for Nuclear Auxiliary Power (SNAP) Reactor — SNAP 
Program," as given in the Standard Distribution for Classified 
Scientific and Technical Reports, M-3679, and according to the 
NASA Supplementary Report Distribution List at the end of this 
report. 
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ABSTRACT 

A digital computer code was developed to simulate the time- 
dependent behavior of the 5-kwe Reactor Thermoelectric System. 
The code was used to determine lifetime sensitivity coefficients 
for a number of system design parameters, such as thermoelec- 
tric module efficiency and degradation rate, radiator absorptivity 
and emissivity, fuel element barrier defect constant, beginning- 
of-life reactivity, etc. A probability distribution (mean and 
standard deviation) was estimated for each of these design param- 
eters. Then, error analysis was used to obtain a probability 
distribution for the system lifetime (mean = 7.7 years, standard 
deviation =1.1 years). From this, the probability that the system 
will achieve the design goal of 5 years lifetime is 0.993. This 
value represents an estimate of the degradation reliability of the 
system. 
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I. INTRODUCTION 


The objectives of the studies described here were threefold: 

1) Develop a model to predict the performance of the 5-kwe Reactor 
Thermoelectric (TE) System over its operating lifetime, from the 
start of full power operation until, through component degradation, the 
system is no longer capable of producing the required power level of 

5 kwe. Use the model to determine the lifetime behavior of the system 
as it is nominally expected to perform in the flight environment. 

2) Use the system model to determine lifetime sensitivity coefficients 
[(3L/L»)/(3X/X)] for each of the key designparameters affecting system 
performance. These coefficients allow the various component per- 
formance parameters to be ranked in the order of their influence on 
system lifetime, and thus provide a basis for allocating design margins 
and for allocating development effort toward those items which are 
most critical to mission success. 

3) Determine a numerical estimate of the degradation reliability of the 
system. A probability distribution is estimated for each parameter 
for which a lifetime sensitivity coefficient has been determined. Using 
the law of propagation of errors, the mean system lifetime is deter- 
mined, as well as the standard deviation in system lifetime. From 
these data, the probability that the system will fail in the degradation 
mode (i. e. , will not, due to component degradation, achieve its design 
lifetime of 5 years) is calculated. One minus this failure probability 
is the system reliability, considering only the degradation modes of 
failure. The overall system reliability, of course, must include the 
catastrophic failure modes as well. Catastrophic failure is not treated 
in this report. 
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II. DESCRIPTION OF MODEL 


A. SUMMARY DESCRIPTION OF SYSTEM AND MODEL 

The 5-kwe Reactor Thermoelectric System is designed to provide electric 
power in a space environment for 5 years. Thermal power is generated by the 
reactor, and transferred to the primary loop NaK coolant. Thermoelectric 
modules convert this thermal power to electrical power. The electrical power 
is used to drive the electromagnetic (EM) NaK pumps for both the primary and 
secondary loops, and to provide 5 kwe at the spacecraft mating interface. The 
secondary NaK loop and the radiator provide heat rejection capability for the 
system. A schematic of the system is shown in Figure 1. 

The system model was developed from the digital computer models used to 
analyze each system component, and is subdivided into a power conversion sys- 
tem (PCS) model and a reactor model. The overall model is schematically 
illustrated in Figure 2. During nominal operation, the temperatures throughout 
the system increase, to compensate for degradation (particularly of the thermo- 
electric modules), and thus produce a constant 5 kwe of power. The end of sys- 
tem life occurs when 5 kwe canTro longer be produced. In the pres-ent model, 
this may happen in one of two ways. Either the reactivity of the reactor becomes 
zero, or module degradation becomes so large that the required power cannot be 
produced, no matter how high the system temperatures are elevated. 

B. POWER CONVERSION SYSTEM MODEL 

The PCS model used in this study, called SYSTEM, is described in detail in 
Reference 1. The components modeled are: (1) the TE power modules, (2) the 

TE pump modules, (3) the radiator, (4) the primary and secondary NaK pumps, 
and (5) the primary and secondary loop piping. The radiator model is multinode, 
whereas all other components are single node. The model utilizes temperature- 
dependent power and pump TE module degradation of both efficiency and power 
per module. 

The SYSTEM code starts by performing a system heat balance to determine 
beginning-of-life (BOL) temperatures around both primary and secondary loops. 
Also determined are the primary and secondary flow rates, and the reactor 
thermal power required. Then, a time step is taken, time- and temperature- 
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INPUT 


DATA CHARACTERIZING ALL COMPONENTS OF 
SYSTEM AND THEIR DEGRADATION WITH TIME 


PCS 

MODEL 


PERFORMS SYSTEM HEAT BALANCE - 
DETERMINES SYSTEM TEMPERATURES AND 
REACTOR POWER AS A FUNCTION OF TIME 
OVER LIFE 


1 

REACTOR 
MODEL 


DETERMINES REACTIVITY AS A FUNCTION 
OF TIME 


OUTPUT 


TEMPERATURES, POWER, AND REACTIVITY 
vs TIME. END-OF-LIFE OCCURS WHEN 
SYSTEM CAN NO LONGER PRODUCE 5 kwe 


6532-5002 


Figure 2. 


Schematic of System Model 
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dependent degradation is computed, and another heat balance is performed. 

This procedure is repeated until a given time has elapsed, or until the system 
can no longer produce the required power of 5 kwe. Thus, a time-dependent 
history of system temperatures, reactor power, and coolant flow rates is ob- 
tained. This part of the model assumes that the reactor can supply any thermal 
power level required, at any reactor outlet temperature required by the system. 
The "system" end of life (EOT) occurs when the module degradation is so large 
that no value of reactor power or outlet temperature will result in 5 kwe output. 
Figure 3 shows, for the nominal system, reactor outlet temperature and reactor 
AT as a function of time. The system EOL, shown as a dashed line, occurs at 
8.8 years. 

C. REACTOR MODEL 

The reactor model used in these studies is based on the parametric reactor 

( 2 ) 

analysis code ZIP. It models the hot fuel element and the average fuel ele- 
ment in the reactor in 1 1 axial nodes, and determines time-dependent fuel tem- 
peratures, fuel swelling, and hydrogen leakage from the fuel. These are based 
on the time-dependent reactor power and inlet and outlet temperatures calculated 
previously by the SYSTEM subroutine. The reactor model then does a reactivity 
lifetime calculation, in which all reactivity losses (due to such factors as hydro- 
gen leakage, fission product buildup, uranium burnout, xenon-135, etc.) are sub- 
tracted from the initial BOL excess reactivity. In this manner, the reactivity 
EOL is determined. For the nominal case, the reactivity EOL occurs at 7.7 years 
(see Figure 3). As will be shown later, either type of EOL may occur first. 

A source program listing of the reactor subroutine and the control program 
is given in the appendix. A listing of the SYSTEM subroutine is given in Refer- 
ence 1. 
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III. ANALYSES PERFORMED 

A. LIFETIME BEHAVIOR OF NOMINAL SYSTEM 

The overall system model was first used to determine the behavior of the 
base-case reactor TE system when operated at 5 kwe, Some of the key system 
parameters for this case are given in Table 1. Note that sufficient margin is 
available to provide a lifetime of 7.7 years, Z.7 years above the required 5. At 
this point, reactivity is zero, and the reactor outlet temperature is 1248° F. 


TABLE 1 

BASE-CASE CHARACTERISTICS 


Parameter 


Value 

at 


BOL 

5 years | 

7.7 years 

Reactor 




Outlet Temperature (°F) 

1122 

1170 

1248 

At ( ° f) 

83.3 

89.4 

99.9 

Power (kwt) 

90.2 

96.2 

106.8 

Thermoelectric Module 




Hot Cladding Temperature (°F) 

1069 

1113 

1188 

Cold Cladding Temperature ( ° F) 

494 

510 

536 

Power Module Degradation 

0.000 

0.051 

0.1061 

Power Module Efficiency (%) 

6.65 

6.24 

5.63 

Hot Reactivity ($) 

2.27 

1.33 

0.00 

Hydrogen Loss ($) 

0.00 

0.47 

0.93 


The first parametric study undertaken was to investigate the lifetime be- 
havior of the nominal system when operated at power levels other than the re- 
quired 5.0 kwe. The results of this study are illustrated in Figure 4, which 
shows time -dependent reactor outlet temperature for power levels between 
4.8 and 5.3 kwe. Also shown in this figure is a dashed line, representing the 
set of points at which reactivity is zero. The same information is shown in 
Figure 5, which is a cross plot of the data of Figure 4. Note that, for powers 
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Figure 5. System Lifetime vs Power 


TABLE 2 

THERMOELECTRIC MODULE SENSITIVITY COEFFICIENTS 


Variable 

Base 

Value 

Sensitivity 

Coefficient 

COL/D/OX/X)] 

Power Modules 



Power per Module (Normalized) 

1.0 

+3.4 

Module Efficiency (Normalized) 

1.0 

+4.8 

Degradation Rate (fraction/year) 

0.010 

-0.30 

Pump Modules 



Power per Module (Normalized) 

1.0 

-0.54 

Module Efficiency (Normalized) 

1.0 

+ 0.63 

Degradation Rate (fraction/year) 

0.014 

-0.001 


Value 
Resulting 
in 5-year 
Lifetime 


0.920 

0.930 

0.019 

1.68 

0.62 

Large 


TABLE 3 

RADIATOR SENSITIVITY COEFFICIENTS 


Variable 

Base 

Value 

Sensitivity 

Coefficient 

[OL/L)/(3X/X)J 

Value 
Resulting 
in 5 year 
Lifetime 

Emissivity 

0.91 

+ 5.3 

0.85 

Solar Absorptivity 

0.40 

-0.27 

0.93 

Fin- Tube Thermal Bond 
Resistance ( ° F) 

2 

-0.039 

17.2 

Radiator Dimensional 
Tolerance Factor 
(Normalized) 

1.0 

+ 3.2 

0.87 
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below 5.15 kwe, reactivity EOL is the limiting parameter. Above 5.15 kwe, 
system EOL is limiting. From this study it was determined that: 

1) System lifetime is sensitive to electric power demand at a rate of 
1.3 years /0.1 kwe 

2) The nominal lifetime margin of 2.7 years would be used up, if the 
system were operated at a power of 5.19 kwe. 

B. DETERMINATION OF LIFETIME SENSITIVITY COEFFICIENTS 

A number of key design variables were selected for sensitivity coefficient 
determination, and the system model was run to obtain each coefficient. In 
many cases, two off-nominal values were run for a given variable, where it was 
suspected that the lifetime-vs-variable curve was nonlinear. In these cases, one 
off-nominal point was selected on each side of the nominal point. 

The results of the lifetime sensitivity coefficient study are shown in Tables 2 
through 5, and in Figure 5. 

Table 2 lists the variables associated with the thermoelectric power and 
pump modules that were selected for coefficient determination. These are the 
electrical power produced by a module for a given hot cladding and cold cladding 
temperature, the module efficiency for a given hot and cold cladding temperature, 
and the module degradation rate per year at an average temperature of 1085° F. 
Note that the power module parameters are much more critical than the pump 
module parameters. For example, reducing the power module efficiency by 7% 
(while holding all other parameters at their nominal values) results in reducing 
the nominal 7.7-year lifetime to the required 5.0-year value. 

Figure 6 shows the strong influence of power module degradation rate on 
reactor outlet temperature. Virtually all of the PCS performance degradation 
over lifetime results from power module degradation (TE module degradation, 
in the model used here, applies both to power per module and module efficiency). 

Table 3 lists the key radiator parameters and their sensitivity coefficients. 
Also given is the value of each parameter which would result in a 5.0-year life- 
time (while holding all other parameters at their nominal values. From this 
table, it may be seen that the emissivity of the radiator surface is a critical 
parameter. 
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TABLE 4 


PUMP ELECTROMAGNETIC AND LOOP 
HYDRAULIC COEFFICIENTS 


Parameter 

Base 

Value 

OL/L)/OX/X) 

Primary Pump Throat 
Magnet Field Strength (G) 

2380 

+0.008 

Secondary Pump Throat 
Magnet Field Strength (G) 

2380 

+ 0.005 

Electrical Resistance of 
Pump Bus and Braze 
Joint (Normalized) 

1.0 

+0.013 

Primary Loop Hydraulic 
Resistance 

0.03576 

j 

-0.01 

Secondary Loop Hydraulic 
Resistance 

0.1519 

-0.004 


TABLE 5 

REACTOR SENSITIVITY COEFFICIENTS 


Variable 

Base 

Value 

Sensitivity 

Coefficient 

COL/D/OX/X)] 

Value 
Resulting 
in 5-year 
Lifetime 

Barrier Defect Constant, 
Ad (Fractional area of 
defects) 

0.0015 

-0.16 

0.0064 

Fuel -Cladding Gap at 
BOL (mils) 

12 

-0.034 

Large 

Hot BOL Reactivity ($) 

2.27 

2.28 

0.94 

Reactivity Loss Rate, 
Excluding Hydrogen Loss 
($/year) 

0.17 

-3.19 

1 

0.36 
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REACTOR OUTLET TEMPERATURE (°F) 



6532-5006 

Figure 6. Reactor Outlet Temperature vs Time, With Power 
Module Degradation Rate (%/year at 1085° F) as Parameter 
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Figure 7. Probability Tree for Pump 
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Table 4 lists the parameters selected for sensitivity coefficient analysis in 
the pump electromagnetics and loop hydraulics areas. Since the coefficients for 
all these variables are small, no column was entered showing the value resulting 
in a 5-year life. 

Table 5 lists the key reactor variables selected for coefficient determination. 
From this table, it may be seen that the BOL fuel- cladding gap has a weak in- 
fluence on lifetime, but that BOL reactivity and the nuclear reactivity loss rate 
are important. Here, the nuclear reactivity loss rate includes uranium burnout 
and fission product buildup. 

C. DEGRADATION RELIABILITY OF SYSTEM 

In general, reliability prediction is accomplished by: (1) construction of an 

analytical model to represent the system, (2) determination of the necessary 
physical constants, failure rates, and probability distributions required for in- 
put to the model, and (3) performing a calculation with the model to find the 
probability that the system will not fail prior to achieving its design objective 
(in this case, producing 5 kwe for 5 years). This overall success probability is, 
by definition, the predicted reliability of the system. 

The standard approach to the problem of reliability prediction is the use of 
a probability tree to characterize all failure mechanisms for each individual 
component of the system. The separate trees for each component are then inte- 
grated into an overall system tree, to obtain a system reliability prediction. A 
quantitative success probability estimate is obtained for each block at the base 
level of the tree, and these probabilities are combined appropriately, through 
and-gates and or-gates, to determine the success probability at each higher 
level. An illustrative example of a probability tree for the pump is shown in Fig- 
ure 7. The probability tree is an adequate model to represent the catastrophic 
failure mechanisms whereby the system fails suddenly. 

For degradation modes of system failure, in which the loss of performance 
is gradual, and may or may not cause failure, the probability tree is not suffi- 
cient. This is particularly true for interactive modes of performance degrada- 
tion; for example, degradation of the radiator emissivity coating results in in- 
creased system temperatures, which then results in increased degradation of 
power module efficiency. To properly account for degr adational failure modes, 
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TABLE 6 

COMBINATION OF UNCERTAINTIES 


Variable 

Base 

Value 

^Variable 

ffLife 

(years) 

Barrier Defect Constant, Aj} 

0.0015 

0.00037 

0.30 

Fuel-Cladding Gap at BOL (mils) 

12 

1.75 

0.07 

Hot Reactivity at BOL ($) 

2.27 

0.25 

0.37 

Reactivity Loss Rate — Nuclear ($/year) 

0.17 

0.06 

0.74 

Radiator Emissivity 

0.919 

0.011 

0.46 

Emissivity Degradation Rate (fraction/year) 

0 

0.001 

0.14 

Solar Absorptivity 

0.3 

0.05 

0.26 

Absorptivity Degradation Rate (fraction/year) 

-0.04 

0.02 

0.12 

Fin-Tube Bond Thermal Resistance (°F) 

2 

0.1 

0.02 

Bond Degradation Rate (°F/year) 

0.05 

0.02 

0.02 

Radiator Dimensional Tolerance 

1.0 

0.0024 

0.06 

Pump Magnet Field (G) 




Primary 

2380 

15 

0.004 

Secondary 

2380 

15 

0.001 

Magnet Degradation (fraction/year) 




Primary 

0.03 

0.01 

0.007 

Secondary 

0.03 

0.01 

0.015 

Pump Bus and Braze Joint Electrical 
Resistance 

1.0 

0.04 

0.050 

Resistance Degradation Rate (fraction/year) 
Hydraulic Resistance 

0.005 

0.0017 

0.005 

Primary Loop 

0.0358 

0.0021 

0.040 

Secondary Loop 

0.1519 

0.0076 

0.003 

Root-Sum-Square Combination of O t 

Life 

's 


1.1 year 
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the system model described in this report is used. Specifically, a probability 
distribution is estimated for each of the key system performance parameters 
affecting system lifetime. For each such parameter, a lifetime sensitivity coef- 
ficient is calculated, using the system performance model, as previously dis- 

(3) 

cussed. Then, the law of propagation of errors is used to obtain an estimate 
of the probability distribution associated with the system lifetime. From the 
mean and standard deviation of this distribution, the probability of not achieving 
the design lifetime due to excessive degradation may be determined. Degradation 
reliability is simply one minus the degradation failure probability. 

Table 6 lists each of the key parameters, gives the expected mean value of 
each, and the estimated one- standard-deviation uncertainty in each. The ap- 
propriate sensitivity coefficients are then used to obtain the final column of 
Table 6, the lifetime uncertainty due to the expected uncertainty in the given 
parameter. Finally, root- sum- square combination of the lifetime uncertainties 
results in a standard deviation in the overall lifetime of 1.1 years. In conjunc- 
tion with the mean lifetime estimate of 7.7 years, this implies a system degrada- 
tion reliability of 0.993. 

Note that the TE module performance parameters are not included in Table 6, 
and thus do not contribute to the degradation reliability value given. The justifi- 
cation for this omission is that the modules are still under development, and the 
exact level of their performance and the uncertainty in this performance, after 
completion of the development program, is virtually impossible to estimate at 
this time. Therefore, it was assumed that the development program will, as a 
minimum, achieve its current performance goals, which would result in module 
performance as good as, or better than, that used in the system model. Thus, 
the system degradation reliability quoted will be achieved or exceeded, if the TE 
module development program is successful. 
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IV. CONCLUSIONS 


The principal conclusion of these studies is that the current reference 5-kwe 
Reactor Thermoelectric System is estimated to have a mean lifetime of. 7.7 ±1.1 
years, which implies a degradation reliability of 0.993 and therefore a 0.7% 
probability of not achieving the design lifetime of 5 years. Although this esti- 
mate is approximate, its accuracy appears adequate for the current conceptual 
design phase. Further, the degree of overall design margin implied by this 
estimate appears to be adequate, and no design changes are recommended at 
this time. 
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I 


- 1 - 


C0NTRL 


1 00* 0VR, SYSTEM 
1 1 O*0VR,REX 
1 20*LIB, INTRP2 
1 30* LIB, BNDHZR 
1 40 INTEGER NCAC10) 

150 REAL EP0WERC2 1 ) 

160 REAL T0UTLTC21 >,PREACTC21 >,C0REDTC21 ) 

170 REAL TH0TCLC21 >,TC0LCLC2i > 

180 REAL EFFP0MC2 I >, DGRP0MC2 1 > 

190 REAL EFFPMMC21 >, DGRPMMC21 > 

200 REAL TAVRADC21 ),PRADTRC2J >,DTRADTC21 ) 

210 REAL T I NLET (21)/ RH0 C 2 1 ) 

220& , H3LC2I >, SNAMEOO) 

230 REAL XDATAC 50, 10) 

240 REAL DUMMY C 25 ) 

25044, TH0TCPC 21 ),TC0LCPC2l ) 

26044, T YME ( 2 1 ),V0LTAGC2J ) 

2708c j BASE (35) 

280 C0MM0N DUMDUM ( 933 >, I 1 1 DUM ( 4 ) 

290 DATA SNAME/"Sr , ,"S2","S3”,"S4","S5","S6", 

30044 "S7", "58", "S9", "S 1 0", "S 1 1 'L "SI 2 "* "S 1 3", 

31 044 "SI 4","S15",”S16","S1 7", "SI 3", "SI 9","S2 0" 

32044 , ”S2 1 ”, "S22", "S23", "S24”, "S25", "S26", ”S27", "S28", "S29", "S30" 

330& / 

340 80 CONTINUE 

350C BEGIN SUPER-L00P, I . E . , NEW FILE CASE. 

360 PRINT, t >"NO. CASES, CASE ID NUMBERS", » * 

370 READC50, >, NCASES, CNCACJ),J * 1, NCASES) 

380 D0 82 J = 1 , NCASES 

390 JG = NCAC J > 

400C TRANSFERS THE DATA FR0M FILE INTO PROGRAM 
410 CALL OPENF C 1 , SNAME C JG > ) 

420 READC 1 , ),NVAR, CXDATAC I, J), I 1,NVAR) 

430 CALL CLOSEF C 1 ) 

440 82 CONTINUE 

45 OC GIVES NAMES T0 FILE DATA. VARIABLES "P0MPPM" THROUGH "RDIMTL" 

46 OC ARE EXPLAINED IN "SYSTEM" ROUTINE. 


470 

D0 83 M 


1, NCASES 

480 

P0MPPM 

3 

XDATAC 1 ,M) 

490 

POMEFF 

= 

XDATA( 2, M ) 

500 

POMDGB 

= 

XDATA( 3, M ) 

51 0 

PMMPPM 

3 

XDATAC 4, M > 

520 

PMMEFF 

3 

XDATA( 5, M ) 

530 

PMMDGB 

3 

XDATAC 6, M) 

540 

PTFSIN 

3 

XDATAC 7, M) 

550 

PTFSDG 

3 

XDATAC 8, M) 

560 

STFSIN 

* 

XDATAC 9, M) 

570 

STFSDG 

3 

XDATAC 1 0,M> 

580 

PBUSIN 

* 

XDATAC 1 1,M) 

590 

PBUSDG 

3 

XDATAC 12, M) 
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2 


CONTRL CONTINUED 


600 PPIPFH * XDATAC 1 3, M) 

610 SPIPFH ■ XDATAC I 4,M) 

620 REMMIN * XDATAC 1 5, M) 

630 REMMDG = XDATAC 1 6, M) 

640 RABSIN * XDATAC 1 7/ M) 

650 RABSDG * XDATAC 1 8, M) 

660 RBNDJN * XDATAC 1 9, M ) 

670 RBNDDG = XDATAC 20, M) 

680 RDJMTL * XDATAC 2 1,M) 

690 NSTEPS » XDATAC 22, M) 

700C MAX. NO. 0F TIME STEPS CEACH STEP*l/2 YEAR) 

710 EP0WB = XDATAC 23, M ) 

72 0C ELECTRICAL POWER CCONSTANT IN TIME) 

730C •’SYSTEM" WILL GIVE P0WER AS A TIME-DEPENDENT VARIABLE. 

740 KSYS “XDATAC 31 ,M) 

750 NW 0HD“XDATAC 35, M ) 

76 OC KSYS=0/1 — QUANTITIES ARE TIME DEPENDENT/CONSTANT 
770C NW0RD IS -2/0/1 “EXECUTES SYSTEM AND SAVES OUTPUT F0R NEXT 
780C ROUTINE IN A FILE/EXECUTES SYSTEM /DOES N0T EXECUTE SYSTEM 
790C BUT READS THE FILECFR0M PREVIOUS CASE) FOR THE NEXT ROUTINE. 

800 DO 31 N=24, 3 0 

810 NZ=N-23 

820 DUMMYCNZ) = XDATAC N,M) 

830C EXPLAINED AT THE BEGINNING OF "REX" 

840 81 CONTINUE 

850 IF CKSYS «EQ • 1 • OR. NW 0RD. EQ. 1 ) GOTO 30 

860 LCNTRL « 2 

870 CALL LI NKC 4, "OSYSTE” ) 

880 CALL SYSTEM CPOMPPM, 

8904 P0MEFF , POMDGB, EP0WB, 

9004 PMMPPM, PMMEFF, PMMDGB, 

91 04 PTFSIN.PTFSDG, STFSIN, STFSDG, PBUS I N, PBUSDG, 

9204 PPIPFH, SPIPFH, 

9304 REMMIN, REMMDG, RABSIN, RABSDG, RBNDIN, RBNDDG, RDIMTL, 

9404 NSTEPS, LCNTRL, KCNTRL, 

9504 EPOWER, VOLTAG, F INTIM, 

9604 TOUTLT, PREACT, C0REDT, 

9704 THOTCL, TC0LCL, 

9804TH0TCP, TCOLCP, 

9904 EFFPOM, DGRP0M, EFFPMM, DGRPMM, 

10004 TAVRAD, PRADTR, DTRADT ) 

10I0C EXECUTES "SYSTEM” AND PROVIDES TIN, TOUT, COREDT, POWER AS 
1 02 OC FUNCTIONS OF TIME. 

1030 PRINT 99, FINTIM 

1040 99 F0RMATC”F INTIM * ",F7.4) 

1050 PRINT, 

1060 I F C NW ORD • NE • -2 ) G OT 0 30 
1070 WRITEC2, 3333), NSTEPS 
1080 LIFE-NSTEPS+1 

1090 WRITEC2, 500), CT0UTLTC I >, COREDT C I ), PREACTC I >, EP0WERC I ), I » 1 ,LIFE ) 
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CONTRL CONTINUED 


1100 WR I TEC 2, 600 ) , F I NT IM 
1110 600 F 0RMATCF 11.4) 

1120 500 F0RMATCC4F1 1 .4)) 

1 130 3333 F 0RMAT (16) 

1140 CALL CL0SEFC2) 

1150C SAVES "SYSTEM" 0UTPUT F0R FUTURE USE <NW0RD*-2 HERE). 

1160 30 CONTINUE 

1170 IF (NWORD.NE. 1 )G0T0 700 

1180 CALL 0PENF ( 3, "MORDAT" ) 

1190C READS THE FILE WITH "SYSTEM" 0UTPUT FR0M SOME PREVIOUS RUN. 
1200 READC 3s 3333* ERR=66 ) , NSTEPS 
121 0 LIFE=NSTEPS+ I 

1220 READ <3,500, ERR=6 7 ) , C T 0UTLT C I ) , C OREDT C I ) , PREACT < I ) , EP0WER Cl), 
1230& 1=1, LIFE) 

1240 READC 3, 600, ERR=65 ),FINTIM 
1250 G0 TO 68 

1260 66 PRINT, »2," ERROR IN NSTEPS READ", t 2 
1270 CALL EXIT 

1280 67 PRINT, t 2, " ERROR IN ARRAY READ", ♦ 2 
1290 CALL EXIT 

1300 65 PRINT, t2, "ERROR IN READING FINTIM",t2 

1310 CALL EXIT 

1320 63 CALL CL0SEFC3) 

1330 700 CONTINUE 
1340 DO 85 JTIME= 1,21 

1350 85 TYME ( JT IME ) =• 5+FLOAT ( JTIME- 1 ) 

1360 LIFE= 1 +NSTEPS 

1370 IF CKSYS.EQ.OIG0T0 35 

1380C IF THE OPTION TO AVOID "SYSTEM” HAS BEEN USED, I.E., IF 
1390C QUANTITIES ARE CONSTANT IN TIME, TOUT, COREDT, REACTOR POWER 
1400C ARE READ IN FROM FILE WHICH SHOULD PROVIDE THEM. 

1410 DO 10 1=1, LIFE 

1420 TOUTLTC I ) =XDATA< 32, M > 

1430 C0REDT Cl) =XDATAC 33, M ) 

1440 PREACT C I > -XDATAC 34, M > 

1450 10 EPOWERCI )=EP0WB 

1460 35 CONTINUE 

1470 DO 84 J = 1 , L IFE 

1480 TINLETC J) = TOUTLTC J ) -COREDT < J ) 

1490 84 CONTINUE 

1500 2300 FORMAT*. " SYSTEM E 0 L IS ”,F10.2) 

1510 CALL L I NK C 5 , " OREX" ) 

1520C EXECUTES THE MAIN PART Tnc. PROGRAM. FINDS REACT 0R 
1530C CHARACTERISTICS AS A FUNCTION OF TIME. 

1540 CALL REXCT INLET, COREDT, PREACT, NSTEPS, RH0, T I MEND, H3L, DUMMY) 
1550 ZER0=O. 

1560 PRINT 2200, ZERO, TIMEND 

1570 2200 FORMATC " TIME (WHEN REACTIVITY IS", FI 0.1," ) * ”,F10.2) 

1580 VALUE= 1 200 
1590 NL0C=O 
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C0NTRL CONTINUED 


1600 TTEMP-FUNCT1 (T0UTLT, TYME, 21 , VALUE,NL0C, 2) 

1610 PRINT 21 00, VALUE, TTEMP 

1620 2100 F0RMATC " TIME (WHEN T0UTLT JS ",F10.1,” ) = ", 

1630&F10.2) 

1640 IF (KSYS .EQ • OJPRINT 2300,FINTI.M 
1650 PR I NT , 

1660 PRINT 2000, 

1670 2000 F0RMAT(2X, "TIME", 2X, "KWE", 3X, "TOUT", 3X, "TIN”, 2X, "DEL T", 

16804 3X, "KWT", 3X, "RHO", 3X, ”$ HL") 

1690 PRINT 2001, (TYME(L),EP0WER(L),TOUTLT(L),TINLET(L), COREDT(L), 

1700* PREACT (L ), RHO (L) , H3L < L ) , L = 1,LIFE) 

1710 2001 F0RMAT((F5.1,F7.2,2I6,2F7.2,2F6.2)) 

1720 PRINT, 

1730 IFCKSY5.EQ. 1 .OR.Ntf ORD.EQ. 1 )G0T0 40 

1740 PRINT, t, "TIME THCL TCCL THCP TCCP PM EFF PM DEG PU EFF PU DEG" 
1750 PRINT 2002, CTYMEC I ),TH0TCL( I ), TC0LCLC I ), TH0TCP< I ),TC0LCP( I ), 

1 76O&EFFP0MC I ) , DGRP0M( I ), EFFPMM( I ) , DGRPMM < I ), 1 = 1 ,LIFE) 

1770 2002 FORMAT < (F4. 1 , 4I5,4F8*4) ) 

1780 PRINT, t, "TIME VOLTAGE TAVRAD DTRADT P RAD" 

1790 PRINT 2400, < TYME ( I ) , V 0LTAG C I ) , TAVRAD ( I ) , DTRADT ( I ) , PRADTR ( I ) , 
1800*1=1, LIFE) 

1810 2400 FORMAT < 

18204CF4. J,3X,F6.1,3X,F6.1,3X,F6.1,3X,F6.2)) 

1830 40 CONTINUE 

1840 PRINT, ♦ , "S-F ILE DATA” 

1850C PRINTS OUT THE FILE THAT WAS USED. 

1860 PRINT 2010, POMPPM, POMEFF, POMDGB, PMMPPM, PMMEFF, PMMDGB, 

1 8704PTFS IN, PTFSDG, STFSIN, STFSDG, PBUSIN, PBUSDG, PPIPFH, SPIPFH, 

1 8804REMMIN, REMMDG, RABSIN, RABSDG, RBNDIN, RBNDDG, RDIMTL, EPOWB, 

1890* < DUMMY < I ) , I * 1 , 7 ), NSTEPS, KSYS, NWORD 

1900 2010 FORMAT (/"POMPPM ”,F8.3," POMEFF ",F8.3, 

19104" POMDGB ”,F8.3," PMMPPM ",F8. 3, /"PMMEFF ”,F8.3, 

19204" PMMDGB ",F8.3,” PTFSIN ",F8.I," PTFSDG ",F8.3,/ 

I 9304"STFS IN ",F8. 1 ,” STFSDG ",F8.3," PBUSIN ",F8.3, 

1 9404 " PBUSDG ",F8. 3, /"PPIPFH ”,F8.5," SPIPFH ",F8.5, 

19504” REMMIN ”,F8.3,” REMMDG ",F8. 3, /"RABSIN ",F8.3, 

I960*" RABSDG ",F8.3," RBNDIN ",F8.1," RBNDDG ",F8.3,/ 

1 9704"RDIMTL ",F8.3," EPOWB ",F8.3,” GAP ",F8.4, 

19804" B0LRH0 ”, F8. 3, /"ACOEFF ”,E9.3," AEXP ",F8.I, 

19904" KPT ", 18,” BUZ ”,F8.4, /"AMO ",F8.5, 

20004" NSTEPS ", 18," KSYS ”, 18," NWORD ", 18) 

2010 BASTIM=8 .24 

2020 IF (NCA(M) *£Q.2O)G0TO 366 

2030 CALL OPENF ( I , "BAS DAT" ) 

2040 READd, ),BASTIM 
2050 CALL CLOSEF(l) 

206 OC BASTIM IS THE LIFETIME OF THE BASE CASE AND IS USED FOR 
2070C THE SENSITIVITY COEFFICIENTS CALCULATIONS (IF A FILE 
208 OC PARAMETER HAS BEEN CHANGED IT WILL PROBABLY CAUSE A CHANGE IN 
2090C LIFE TIME! HENCE THE SENS. COEF. -DELTA LIFE/DELTA PARAMETER) 
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C0NTRL CONTINUED 


2100 CALL 0PENF ( 1 , "S20" ) 

2i 10 READC 1, >,NVAR, (BASECI), I*1,NVAR> 

2120 CALL CL0SEF ( 1 ) 

2130 DC 333 1*1,34 

2140 DELTA*XDATAC I,M)-BASE< I > 

2150 JF CABS C DELTA) • GT. 0. 000001 )G0TO 334 

2160 G0T0 333 

2170 334 JACK* I 

2180 GCT0 335 

2190 333 CONTINUE 

2200 335 DELT I M*T IMEND-BAST IM 

2210 C0EFF=DELT I M/DELTA 

2220 PRINT 336, JACK# C0EFF 

2230 336 FORMAT (/” CHANGING THE MEMBER N0. ",I3, 

2240**’ 0F THE XDATA ARRAY” /"GIVES ”,F12.4,/" AS THE 
2250* SENSITIVITY COEFFICIENT (DEL LIFE/DEL X)”) 

2260 IF CNCACM) .NE.2O)G0T0 399 
2270 366 CONTINUE 
2280 WRITEC4, ),TIMEND 
2290 CALL CL0SEF < 4, "BASDAT" ) 

2300 399 CONTINUE 
2310 IF (KSYS»EQ»O)GOT0 83 
2320 PRINT, t2 
2330 PRINT 2016, 

2340 PRINT 2004, CXDATAC I, M) , I =32, 34 ) 

2350 2004 F0RMATC5F14.6) 

2360 2016 FORMAT C" TOUTLT COREDT PREACT") 

2370 83 CONTINUE 

2380 GOTO 80 

2390C END SUPER-LOOP. 

2400 END 
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REX 


30000 SUBROUTINE REX(TINXz DELTXz POWXz NXXz RH0z TIMENDz H3Lz DUMMY) 

30010 REAL TN( 11)/ TCG (ll)zTF(ll)/ TFS ( 1 1 )z TC ( 1 1 ) z 
30020* T I NX (21)* PH I A ( 1 1 ) z 

30030* DELTXC21 >/P0WX(21 >zTMX(2! > z BUX (25)/ BUM < 2 1 )/RH0(21 ) 

30040* / H3LC2 1 ) » DUMMY ( 25 ) / TGS ( 1 1 > 

30050*4 L00G (23 )/MIT(23) 

30060 COMMON PO, ZADz EADz AZRz ALIFz OGPz GOLD/ TFMXz TFMEz CBLz CEL, AHLz 
30070* AKA> AK1 / AK2/WHY/ BSRCH/ TMNF/ BTB0L/ BTE0L/ SADC# TCCE/ SADT/ TCT/ 

30080* SADH/ TCHE/ ZADC/ ZADT/ ZADH/ GTH/ AGL/ CLTH/ CGL/ Z0Pz COLD/ AC 0EF/ 

30090* AEXP/BUZ/P0WJNT/TJN/DELT 

30100* / TCDGPC 1 1 )/TIM(25)z TFUELC 1 1 )/G(l 1 ) z ATFU ( 25 ) / ADCE < 5 ) z 

30110* ADT(5)/ADHE(5)zTBD(5)/AM( 1 0) z TPPKC 25 > / P< 25 ) z 8TNA< 1 1 )z 

30120* BTCDG ( 1 1 >/BDTGP< 1 1 )/BTFUE( 1 1 )zE0TN< 1 1 )/E0TC< I 1 )zE0DT(ll >/ 
301 30* E0TF ( 1 1 ), TNAKC I 1 ),DTGPC 1 1 )zAHZR(25> 

301 40&z H2L(2 1 )/ HYLC25 )z HRCC(25 >/ TSWELLC25 )z G00LC23 ) 

30150* z 0SSWC25 >z BUSVC25 >z T JNZ(25)z DELTZC25 )z P0WZ<25 ) 

30160* zAMl (25 >z AM2 ( 25 >z AM3C25 ) 

30170* zTCLAD(25z 1 1 ) 

30180* zKPTzNJKz JSTPzNSTEPS 

301 90C 

30200 GAP*DUMMY(J) 

3021 OC RADIAL GAP IN MILS 
30220 B0LRH0*DUMMY ( 2 ) 

30230C B 0 L REACTIVITY IN S 
30240 AC0EF =DUMMY( 3 > 

30250C 0FF *AC 0EFF *EXP(AEXP/TFU£L> 

30260 AEXP*DUMMY(4) 

30270 KPT=DUMMY ( 5 > 

30280C -1 /O/l /2»STANDARD/FULL/MIDI /MINI PRINT 
30290 BUZ* DUMMY ( 6 > 

30300C CORE AVERAGE 8URNUP(MA/0> 

30310 AMO*DUMMY( 7 ) 

30320C INITIAL AM(I) IN H.L. EXPRESSION BELOW 
30330C 

30340 CALL 0PENF ( 1 / "AXDAT" ) 

30350 READdz ) ALIFEz ELNOz ALENz AKWKGz DIACLz DIAFUz AZRz BSRCH 
30360C ALIFE x LIFE IN YRS • 

30370C ELNO t # 0F ELEMENTS 

30380C ALENs FUEL ELEMENT LENGTH IN IN. 

30390C AKWKG: KWT/KG-U 

30400C DIACLt UNC0ATED CLADDING ID IN IN. 

3041 OC DIAFU: FUEL MEAT DIA IN IN. 

30420C AZR: INITIAL H/ZR 

3043 OC BSRCH* 0/l*CALC. MAX. INITIAL H/ZR FOR ALL DELTA FUEL/BYPASS 
30440* zCLTHzCCLDzZ0PzWHYzNJK 
30450C CLTH: CLAD THICKNESS IN MILS 
30460C CCLD CLAD DIFF. CONST. 

30470C Z0P* LE.99*AD DEPENDS 0N AMO FR0M S-FILES/GT.99/LE. 1 99* 

30480C USE LINEAR AD CHEz CEzT)/GT. 1 99*USE EXPONENTIAL TYPE AD 
30490C WHY: HYDROGEN WORTH (S/.1NSUB H> 
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HEX C0NTINUED 


30500C NJKs 0/1 -B0NZER-SWENSON/S8DR CORRELATION (NOT USED) 

305 10& * AK 1 * AGL* CGL* GTH* AK2 

30520C AK1 * CGL* GTH* AK2 ARE USED IN THE HYDR0GEN L0SS EQUATION BEL0W* 
30530C H. L. -CONST. *(AM(I )*CCLD*SQRT(P)*EXP(-AKI /T SUB 0/CLTH+ 

30540C AGL*CGL*P*EXP( -AK2/T SUB C)/GTH) 

30550C I IS THE AXIAL TEMP. INDEX 
30560& * ADCE* ADT* ADHE* TBD 

30570C ADCE (COLD END ) * ADT < TUBE ) * ADHE (HOT END) I SEE ZOP ABOVE 
30580C TBD: TIME FOR AD ARRAYS 
30590& *TMX*TIM 

30600C TMX: ARRAY OF OPERATING TIME STEPS 

3061 OC TIM: TIME ARRAY (MORE DETAILED AT BEGINNING THAN TMX) 

30620& * ZADC* SADC* TCCE* ZADT* SADT* TCT* ZADH* SADH* TCHE 

30630C AM(1 ) =ZAD+SAD( 1 -EXP( DELTA TIME/TC) 

30640C THIS EQUATION APPLIES FOR AM(l) AT THE COLD END 

3065 OC (CE VALUES)* AT THE TUBE (T)* AT THE HOT END (HE) AND USES 

3066 OC ZAD*SAD*TC AT CE* HE# T 

30670& * TN* TCG* TF * TF S* TGS* OGP* OPT* AMARG* PHIA 

306 8 OC TN: NAK COOLANT TEMP. 

30690C TCG: CLAD/GLASS INTERFACE TEMP. 

30700C TF: FUEL AVG. TEMP. 

3071 OC TFS : FUEL SURFACE TEMP. 

30720C 0GP: - 1 /O/l -CONSTANT GAP/CALC. GAP/OPTIMIZES GAP WITH 
307 3 OC AMARG (BELOW) AT E 0 L 

30740C OPT: 0/1 -PRINT GAP (ONLY IF OGP.GE. 0) /BYPASS 
30750C AMARG: USED ABOVE IN OGP 

30760C PHIA: AXIAL TEMP. DISTRIBUTION (BOTTOM TO TOP) 

30770& * FPK* UBUKK* SMK* BUM0M1 * BUM0M2* PPBUK 

30780C THE PARAMETERS BELOW ARE CONSTANTS IN THE FOLLOWING 
30790C REACTIVITY LOSS EQUATIONS: 

30800C FISSION PRODUCTS ( I )=FPK*BURNUP( I ) 

3081 OC U BURNOUT ( I ) =UBUKK*BURNUP( I ) 

30820C SAMARIUM( I ) -SMK0*( 1 -EXP(BURNUP( I )*SMK) ) 

30830C PREP0IS0N BU( I ) -PPBUK* (F RAC 1 *EXP( -BURNUP( I ) +BUM0M1 *2 1 9. 05 ) 

30840C +FRAC2*EXP(-BURNUP(I)*BUM0M2*219.05))* WHERE 

30850C 219.05 IS A UNITS CONVERSION FACTOR* INDICES I AND 2 

30860C REFER TO THE TWO POISONS IN THE REACTOR 

30870& * FRAC 1 * FRAC2* TDEFK* XENONK* HRDK, PDEFK* SMKO 

30880C XENON ( I )— XEN0NK*POW£R( I ) 

30890C TEMPERATURE DEF£CT( I )-TDEFK*TEMP( I ) 

30900C POWER DEF . C I ) =PDEFK*POWER( I ) 

3091 OC H DEF. (I )-HRDK*POWER( I > 

30920C THE INDEX I REFERS TO THE TIME STEPS 
30930 CALL CLOSEF(I) 

30940C 

30950 NSTEPS-NXX 

30960 GOLD-GAP*! 000. 

30970 POWINT-O* 

30980 BUZ-BUZ/Cl 10. *8766. *5) 

30982 P0WXX=P0WX(21 ) 
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REX CONTINUED 


30983 DELTXX® DELTXC 2 1 ) 

30984 TINXX=TINX( 2 1 > 

30990 CALL TEMPS CGAP, ALIFE, DIAFU, DIACL, ELN0, ALEN, PH JA, 

31000* TINX* DELTX* POWX/ TMX# POW, BUX, TN, TCGj TF » TFS* TGS# 0PT * AMARG ) 

31 01 0C USES TIN* T0UT T0 FIND AXJAL TEMPERATURE DJ STRIBUTI 0NS 

31020 IFCBSRCH.GT.O) G0 T0 12 

31030 AZZ=AZR 

31 040 TMMMF=TMNF 

31 050 ATFUEL=ATFU< 1 ) 

31060 CALL BNDHZRCBBZR* AZZ* ATFUEL* TMMMF ) 

31070C FINDS H/ZR 
31080 AZR-AZZ 

31090 12 CONTINUE 

31 100C 

31110 CALL BHLOSSCALIFE*ELN0*ALEN*DIACL*AHL0S*DIAFU*AMO> 

31120C HYDROGEN LOSS CALCULATIONS 
31130 CALL BETABN(JSTP) 

31140C BETA BOUNDARY CALCULATIONS 
31150 GGP=GAP*1 000. 

31 J60C 

31170 H2L ( 1 ) =»0 

31180 BUMC 1 > = 0» 0 

31190 TINXC 1 >*TINZC 1 ) 

31200 DELTXC 1 >»DELTZC 1 ) 

31210 P0WXC 1 >=POWZC 1 ) 

31220 DO 45 KK=2*21 

31230 JJ=KK+2 

31240 H2LC KK> =HYL( J J ) 

31250 BUMCKK)=BUXC JJ) 

31260 T I NX C KK.) *T I NZ ( JJ) 

31270 DELTXC K.K) = DELTZCJJ) 

31280 POWX (KK)=*POWZ( JJ) 

31290 45 CONTINUE 

31300 T I NX (21 ) *T I NXX 

31301 DELTXC21 )»DELTXX 

31302 P0WXC21 )=P0WXX 

31310 DO 43 LL=1*21 

31320 H3L(LL>-H2L(LL) 

31330 48 CONTINUE 

31 340C 

31350 CALL XPLODE C T I NX. DELTX* TMX* BUM* POWX* RH0* T IMEND* B0LRH0 

31 360* *FPK*UBUKK* SMK* BUMOM 1 * BUM0M2* PPBUK 

31 370* * FRAC 1 * FRAC2* TDEFK* XENONK* HRDK* PDEFK* SMKO 

31380* ) 

31390C BALANCES REACT IV IT IES* F INDS EOL TIME (WHEN SUM 0F REACTIVITIES 
31400C EQUALS ZERO) 

31410 IF ( KPT «LT. 0 ) GOTO 47 
31420 JMZ* JSTP+ 1 

31430 JMIN=JSTP 

31440 JMAX= JSTP 
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REX CONTINUED 


31450 JF(KPT.LT.I) JMIN-1 

31460 PRINT 31, 

31470 31 F0RMAT<2X, “TIME",2X,"T IN",2X,"DEL T'*, 2X, "P0W", 3X, 

314804 "M A 3X, "TSWEL", 2X, “0Fi* Sr.T**, 2X, u bU SW", <*X, "GAP" ) 

31490 PRINT 32, <TIM< J>, TINZC J) , DELTZ C J >, P0WZC J),BUX< J), 

31 5004 TSWELL( J), 0SSW< J),BUSW( J), G00LC J), J=*JMIN, JMAX) 

31510 32 F0RMAT<I6,F6.1,2F7.1,F8.4,F6.1,3F8.4> 

31520 PRINT 33,TIM<JMZ),BUX<JMZ), 0SSWC JMZ),BUSV< JMZ), G00LC JMZ) 

31530 33 F0RMATCI6,2OX,F8.4,6X,3F8.4) 

31540 PRINT, 

31550 PRINT 85, 

31560 PRINT 86, <TIM< J), AHZRC J), ATFUC J), PC J), HRCC ( J) ,HYL( J>, 

315704 AMI <J), AM2<J),AM3<J), J*JMIN, JMAX) 

31580 PRINT 87,TIM< JMZ ), AHZRC JMZ >, P< JMZ),HYL< JMZ) 

31590 85 F0RMAT<2X,”TIME”,3X,”H/ZR",2X,”AV TF”, 2X, "PRESS”, 

316004 4X,”CC/HR”,2X,”$ H L0SS”, 2X, ”AD1 ”, 5X, ”AD2”, 5X, ”AD3” ) 

31610 86 F0RMATCI6,F8.4, 16, 2F9 . 5, F6 . 2, 3F8 . 4 > 

31620 87 F0RMATC I6,F8.4,6X,F9.5, 9X,F6.2) 

31630 PRINT, 

31640C 

31650 IF(KPT.GT.O) G0 T0 74 

31660 PRINT, 

31670 PRINT 76 

31680 PRINT 77, (BTNAC I ), BTCDG ( I ) , BDTGPC I ) , BTFUEC I), 1*1, 11) 

31690 PRINT, 

31700 PRINT 76 

31710 PRINT 77, (E0TNC I ) , E0TC ( I ) , E0DT Cl), E0TF Cl), 1 = 1, 11) 

31720 PRINT, 

31730 PRINT 76 

31740 PRINT 77, CTNAKC I),TCDGPCI), DTGPC I ) , TFUEL (I), 1 = 1, 11) 

31750 76 F0RMATC5X, "TNAK”, 2X, "CLAD/GL”, 3X, "GAP DT",3X,”AVG FU" ) 

31760 77 F 0RMAT < 4F9.2) 

31770 74 C0NTINUE 
31 780C 

31790 PRINT, 

31800 PRINT 35, P0W, ALIFE, BUX C JMZ ), TIN, DELT, ELN0, GGP, ALEN, DIAFU, 

318104 DI ACL, TFMX, TFME, CBL, AHL, G0LD, AHL05, AZR, Z0P, EAD, P0, VHY, 

318204 BTB0L, TPPKC 1 ),BTE0L,TPPKC JSTP) 

318304 , AK1 , CLTH, AGL, CGL, AK2 

318404 , CCLD 

31850 35 F0RMATC 5X, "P0WER", 2X, "L IFE( YR)", 2X, ”BU (MA/0) ", 5X, 

318604 " I NLET", 2X, "DEL TEMP" /2F 1 0. 1 ,F 1 0.4, 2F 1 0.2/ 

318704 4X, "ELM* TS”, 7X, "GAP", 4X, "LENGTH", 2X, "DIA FUEL", 

318804 2X, "DIA CLAD" /F 1 0.0, 2F1 0.2, 2F 1 0.4/ 

31 8904 4X, "TFMX: 0", 4X, "TF :E0L", 2X, "CC/HR(0>", 2X, ”H/ZR:E0L”, 

319004 6X,”LIFE”/2F10.2,2F10.4,F1 0.2/ 

3 1 9 1 04 1 X, "H L0S S C S > ", 3X, "H/ZR < 0 > ", 3X, "AD C END ) ", 2X, "AD < TUBE ) ", 

319204 2X,"P0<PSIA>"/F1O.2,F1O.3,2F1O.5,F1O.2 

319304 /1X,”W0R<*/NH)",4X, "B0L BT",2X, "B0L TMAX", 4X, "E0L BT E0L TMAX", 

319404 /F 10.3, 4F 1 0.2/4X, "Q CLAD", 3X, "CLAD TH", 2X, "GLASS AG", 
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HEX CONTINUED 


31950* 1 X, "GLASS C0N”,3X,”Q GLASS”/ 

319604 F 1 0.0, F10.1,F10.4,F10.5,F 10.0/ 

31970* 2X, ”CLAD C0N”/F1O.O/) 

31980C 

31990 47 CONTINUE 
32000 LUMP=3+NSTEPS 
32010 PRINT 198, 

32020 D0 101 J= 1 , LUMP 
32030 TIMPAN»TIM( J>/8766. 

32040 101 PRINT 199, T IMPAN, BUX( J) , 0SSV C J ) , BUSW( J ) , G 00L( J) , PC J) 

32050 PRINT 197, 

32060 D0 102 «J“ 1 , LUMP 
32070 TIMPAN»TIM( J>/8766. 

32080 102 PRINT 196, T IMPAN, TSWELL( J), AHZR( J), ATFU( J) 

32090 199 F 0RMAT(F6 • 2, 5F9.4 ) 

32100 196 F0RMAT(F6.2,F9.1,F9.4,F9.1> 

32110 198 F 0RMAT ( /” TIME MA X OFFSET BU SW GAP PRESS”) 

32120 197 FORMAT (/” TIME T SWELL H/ZR AV TF”) 

32130 GAP0* 1 . 

32140 CALL F INDGACG00L, L00G, TIM, MIT, LUMP, GAPTIM, GAP0) 

32150C EXTRAPOLATES T0 FIND THE TIME WHEN THE GAP * GAP0 ( SOME 
32160C SPECIFIED GAP VALUE). 

32170 PRINT 3001, GAP0, GAPTIM 

32180 3001 F ORMAT < /”T I ME (WHEN THE GAP IS ”,F10.1,” ) « ”,F10.2> 

32190 RETURN 
32200 END 
3221 OC 
32220C 
32230C 

32240 SUBROUTINE TEMPS(GAP, ALIFE, DIAFU, DIACL, ELNO, ALEN, PHIA, 

322504 T I NX, DELTX, POWX, TMX, POW, BUX, TN, TCG,TF, TFS, TGS, OPT, AMARG ) 

32260C 

32270 REAL TINXC 2 1 ) , DELTX<2 1 ) , P0WX<2 1 ), TMX (2 1 ), 

322804 PHIA(ll) 

322904 , TNC2 1 ),TCG(21 ), TF (21), 

323004 TFS(21 ),TGS(1 I ),BUX<25) 

3231 OC 

32320 COMMON P0, ZAD, EAD, AZR, ALIF, OGP, GOLD, TFMX, TFME, CBL, CEL, AHL, 
32330* AKA, AK1 , AK2, WHY, BSRCH, TMNF, BTBOL, BTEOL, SADC, TCCE, SADT, TCT, 

32340* SADH, TCHE, ZADC, ZADT, ZADH, GTH, AGL, CLTH, CGL, ZOP, CCLD, ACOEF, 
323504 AEXP, BUZ, POWINT, TIN, BELT 

323604 , TCDGP ( 1 1 ) , TIM(25 ), TFUEL( 1 1 ),G(1 1 ), ATFU(25), ADCE<5 ), 

32370* ADT <5 ) , ADHE(5 ), TBDC5 ), AM( 1 0), TPPKC25 ), P(25 ) , BTNA( 1 1 >, 

323804 BTCDG ( 1 1 ),BDTGP(1 1 ),BTFUE<1 1 ),E0TN(1 1 >,E0TCC1 1 ),E0DTC11 ), 

323904 E0TF ( 1 1 >,TNAK( I 1 >,DTGPC 1 1 ),AHZR(25) 

32400* , H2LC2 1 ), HYL( 25 ), HRCC (25 ), TSWELL(25 ), G00L(23 ) 

324104 , OSSW (25 ), BUSW (25 ) , TINZ(25 ), DELTZ ( 25 ), P0WZ(25> 

32420* ,AMI (25 ) , AM2 (25 ) , AM3 (25 ) 

324304 , TCLAD(25, 1 1 ) 

324404 , KPT, NJK, JSTP, NSTEPS 
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REX CONTINUED 


32450C 

32460 GP=GAP 

32470 AL JF *AL IFE*8766 • 

32480 98 CONTINUE 
32490 JSTEP* 1 

32500 GAP=GP 

32510 GG*GAP 

32520 P0WINT=O.O 

32530 BUX( 1 ) = 0 
32540 D0 17 1=1,11 
32550 17 GCI)=GAP 
32560 99 CONTINUE 
32570C 

32580 TTIM=TIMC JSTEP) 

32590 NLIFE=NSTEPS+1 

32600 NL0C=O 

3261 0 T I N=FUNCT1 CTMX, T I NX* NLIFE, TTIM* NL0C, I ) 

32620 TINZ ( JSTEP )=T I N 

32630 NL0C-O 

32640 DELT=FUNCT1 CTMX, DELTX, NLIFE, TT IM, NL0C, 1 ) 

32650 DELTZ C JSTEP ) =OELT 

32660 NL0C=O 

32670 P0W-FUNCT1 ( TMX, P0WX, NLIFE, TTIM, NL0C, 1 ) 

32680 P0WZC JSTEP )=P0W 

32690C 

32700 JSTEP7* JSTEP+ 1 

32710 P0V I NT=P0W I NT + P0W*C TIMC JSTEP7 ) -T IM( JSTEP) ) 

32720 BUX (JSTEP7) =P0W I NT+BUZ 

32730 AKWFT=P0W*12./CELN0*ALEN) 

32740 TFM=0 • 

32750 TNAKC1)=TJN 

32760 DMP = 0.9752/DIACL 

32770 GMP = GC 1 >*DMP/0.0067 

32780 AKWMP=AKWFT*85.*16./C95.*12. > 

32790 TCDGPC 1 >=TIN+CTCGC 1 J-TNC 1 > )*AKWMP*DMP 

32300 TCLADC JSTEP, 1 >=TCDGPC 1 ) 

32810 DTGPC 1 >=CTFSC 1 I-TGSC 1 ) >*AKWMP*GMP 

32820 DTGL= ( TGS ( 1 >-TCG< 1 > > *AKWMP*DMP 

32330 TFUELC 1 >=TCDGPC 1 J+DTGPC 1 >+CTFC 1 >-TFSCl ) ) *AKWMP+DTGL 
32840 TMNF =TCDGPC 1 ) +DTGPC 1 )+DTGL 

32850 DTM=D£LT/80. 

32860C 

32870 D0 1 0 1=2, 1 1 

32880 TNAK( I > = C TNC I > -TNC I - 1 ) ) +DTM+TNAKC I - 1 ) 

32890 TCDGPC I >=TNAK< I ) + C TCG C I J-TNCI ) )*AKWMP*DMP 
32900 TCLADC JSTEP, I >=TCDGPC I ) 

32910 GMP=GC I >*DMP/0.0067 

32920 DTGPC I > = C TFS C I > -TGS C I ) >*AKWMP*GMP 

32930 DTGL= C TGS C I ) -TCG ( I ) ) *AKtfMP*DMP 

32940 TFUEL C I ) =TCDGP C I > +DTGPC I ) + C TF C I > -TFS C I > ) +AKWMP+DTGL 
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REX CONTINUED 


32950 10 CONTINUE 
32960C 

32970 AT=TFUELC 1 )+TFUEL( 1 1 ) 

32980 KST» 1 

32990 D0 20 1=2, 1 0 

33000 AT =AT+TFUELC I ) *2 . 

33010 TFM=AMAXJ CTFM, TFUELCI)) 

33020 IFCTFUELCII.LT.TFM) GO T0 20 

33030 KST=KST+ 1 

33040 20 CONTINUE 
33050C 

33060 TPPKC JSTEP ) *TFUEL ( KST ) + C TF C KST ) -TFS C KST ) ) *AKWMP 

33070 ATFU < JSTEP ) =AT*0 • 05 

33080 IFC JSTEP. EQ. 1)TFMX=TFM 

33090 IFCJSTEP.NE. 1 )G0 T0 59 

331 OOC 

331 10 DO 59 IJ=1, 1 1 

33120 BTNA ( I J ) =TNAK ( I J ) 

33130 BTCDGC IJ)=TCDGPC IJ> 

33140 BDTGP( I J ) =DTGP< I J ) 

33150 BTFUEC I J)=TFUEL< I J) 

33160 59 CONTINUE 

33170 IFC JSTEP. NE. 1 0)G0 T0 68 

33180C 

33)90 D0 68 I J* 1 , 1 1 

33200 EOTN ( I J ) =TNAK( I J ) 

33210 E0TCC I J)=TCDGPC I J) 

33220 E0DTC I J ) =DTGP( I J) 

33230 E0TF C I J ) *TFUELC I J ) 

33240 68 CONTINUE 
33250C 

33260 CALL GAG C 0PT, JSTEP, GG, GP, BUXC JSTEP7 ), DIAFU,AMARG, PHIA) 

33270C FINDS THE GAP BETWEEN CLAD AND FUEL. 

33280 JSTP*JSTEP- 1 

33290 IFC0GP.LT. O)G0 T0 16 

33300 IF CTIMC JSTEP) .LT.ALIF )G0 T0 99 

33310 IFCGG.NE.GP5G0 T0 98 

33320 16 CONTINUE 

33330 TFME*TFM 

33340 RETURN 

33350 END 

33360C 

33370C 

33380C 

33 3 90C 

33400 SUBROUTINE BHL0SSCAL JFE, ELN0, ALEN, DIACL, AHL0S, DIAFU, AMO) 

3341 OC HYDROGEN LOSS ROUTINE 
33420 DIMENSION TCC10) 

33430 COMMON P0, ZAD, EAD,AZR, ALIF, 0GP, GOLD, TFMX,TFME, CBL, CEL, AHL, 

33440* AKA, AK1 , AK2, WHY, BSRCH, TMNF, BTBOL, BTEOL, SADC, TCCE, SADT, TCT, 
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REX CONTINUED 


33450& SADH^TCHE, ZADC, ZADTj ZADH, GTH, AGL, CLTH, CGL, Z0P., COLD, AC0EF, 

33460& AEXP,BUZ,P0WINT,TIN,DELT 

33470& j TCDGPC 1 1 > , T IM C 25 ) * TFUEL ( 1 1 ),GC J 1 ), ATFU(25>.» ADCE<5>, 

33480& ADT C 5 >.* ADHE C 5 ) * TBDC 5 ) » AMC 10)j TPPKC 25 ># PC 25 ).»BTNA( J I ), 

33490& BTCDGC 1 1 >,BDTGP( I 1 ),BTFUEC 1 1 ),E0TNC 1 1 >,E0TCC 1 1 >,E0DT< 1 I )« 

33500& E0TFC1 1 ) * TNAKC 1 1 >,DTGPC1 I >,AHZRC25> 

33510& j H2L C 2 1 ) , HYLC 25 > , HRCC ( 25 > , TSWELL < 25 > , G00L C 23 ) 

33520& , 0SSW C25 ) .* BUSW C25).TINZC25>* DELTZC 25 ) » P0WZC25 ) 

33530& » AM 1 C 2 5 ) j AM2 C 2 5 ) j AM 3(25) 

33540& j TCLADC 25* 1 1 ) 

33550& *KPT*NJK, JSTP*NSTEPS 

33560C RAYMOND* S PRESSURE EQUATION 

3357 0 F PRES C X* AT ) 3 EXP < - 8 • 8455 + 88 *9801 *X-78. 89t>l *X**2+21 .3731 *X**3> 

33580& *LXPC 1-12. 972 + 9. 770 / *a-2 . *»984*X**2 >*1 .0E04/AT) 

33590 D0 777 JACK=i*10 
33600 777 AMCJACK)=AM0 
33610 JKD=2 

33620 AKA=AGL*CGL/( GTH* . 0254 ) 

33630 J=1 

33640 ZAD= Z0P 

33650 Z 0P=AM Cl) 

33660 EAD=AMC2) 

33670 BKD=CCLD/ C CLTH*0 • 0254 ) 

33680 AHZR Cl) =AZR 

33690 AREA=ALEN*3 .141 59*DIACL*2 .54**2 *0.1 

33700 W0RH=WHY/.O287 

33710 CCELM=22 .4E03 *6 • 3E22*0 .7854*16. 3 9*ALEN*DI AFU**2/( 1 . 8 1 *6 . 023E23 > 

33720 CCELM=CCELM/2. 

33730 D0LCC=W0RH/CCELM 

33740 CC=0.0 

33750 DNH=0.0 

33760 AT=ATFU C 1 ) +45 9.7 

33770 PC 1 ) =F PRES C AZR* AT ) 

33780 P0=PC 1 >*14.696 

33790 IFCZAD.GT. 199. )G0 T0 70 

33800 IFCZAD.GT. 99. >G0 T0 22 

33810 28 C0NTINUE 

33820 ALIF =ALIFE*8766 • 

33830 D0 32 1=1*10 

33840 TCCI >=CTCLADC J* I )+TCLADCJ* 1 + 1 ) >*.5+459.7 

33850 CC=CC+ C AMC I > *BKD*SQRT C PCU > ) *EXPC -AK1 /TCC I ) ) ♦ 

33860& AKA*PC J > *EXPC -AK2/TC C I > > >*AREA*TIMC2) 

33870 32 C0NTINUE 

33880 DNH=CC/CCELM 

33890 CCHR=CC/T IMC 2 ) 

33900 HRCC C 1 ) =CCHR 

33910 CBL=CCHR 

33920 DL0S=DNH*W0RH 

33930 HYLC 1 >=0 

33940 HYLC 2 ) =DL0S 
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REX CONTINUED 


33950 J=2 

33960 87 CONTINUE 

33970 IFCTIMC J+l > . GT . AL IF >TIM< J+ 1 )=ALIF 

33980 AHZRC J>=AHZRCJ-1 J-CC/CCELM 

33990 AZR=AHZRCJ> 

34000 IF ( 0GP.LT. 0)ATFU(J)=ATFU( 1 ) 

34010 AT=ATFUCJ>+459. 7 

34020 PC J)=FPRESCAZR,AT> 

34030 IFCZAD.GT. 199. >G0 T0 70 

34040 IFCZAD.GT. 99. )G0 T0 22 

34050 29 CONTINUE 

34060 TIME=TIMC J+D-TIMCJ) 

34070 81 CC=0. 

34080 D0 57 I = I > 10 

34090 TCC I ) = CTCLADC J, I ) +TCLADC J.» 1 + 1 ) ) *.5+459. 7 

34100 CC=CC+CAMC I ) *BKD*SQRT C PC J ) ) *EXPC -AK1 /TCC I ) )+ 

341 1 0& AKA+P C J ) *EXPC -AK2/TC C I ) ) ) +AREA+TIME 

34120 57 CONTINUE 

34130 DNH=DNH+CC/CCELM 

34140 DL0S»DNH*V0RH 

34150 HYLC J+l >=DL0S 

34160C 

34170 CCHR*CC/TIME 

34180 HRCC < J > =CCHR 

341 90 IF (TIMCJ+1 ) »EQ • AL IF )G0 T0 31 

34200 J=J+1 

34210 G0 T0 87 

34220 70 CONTINUE 

34230 TMM=0.5*CTIMC J+l )+TIMC J) > 

34240 AMC 1 )=SADC*C 1 -EXPC -TMM/TCCE > > + ZADC 

34250 AMC2)=SADT*C 1 -EXPC -TMM/TCT ) I+ZADT 

34260 AMC 1 0 ) =SADH*C 1 -EXPC -TMM/TCHE ) ) +ZADH 

34270 AM 1 C J ) =AM C 1 ) 

34280 AM2CJ)=AMC2> 

34290 AM3C J)=AMC 1 0> 

34300 D0 76 JIL=3,9 

34310 76 AMC JIL)=AMC2) 

34320 IFCJ.EQ.DG0 T0 28 

34330 G0 T0 29 

34340 22 CONTINUE 

34350 IFCTBDC JKD) .GT. TIMCJ+1 ) >G0 T0 13 

34360 JKD= JKD+ 1 

34370 13 TSTP=TBDC JKDI-TBDC JKD- 1 ) 

34380 SC- C ADC E C JKD ) -ADCE C JKD- 1 >)/TSTP 

34390 ST= C ADT C JKD ) -ADT C JKD- 1 ) ) /TSTP 

34400 SH=*CADHEC JKDI-ADHEC JKD-1 )>/TSTP 

34410 24 TMKD*0 »5*C TIMCJJ+TIMCJ+1 ) ) -TBDC JKD- 1 ) 

34420 AMC 1 )=ADCEC JKD-1 J+SC+TMKD 

34430 AMC2)=ADTC JKD-1 >+ST*TMKD 

34440 AMC 1 0 ) =ADHEC JKD- 1 >+SH*TMKD 
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REX CONTINUED 


34450 AM 1 C J > *AM Cl) 

34460 AM2 ( J ) =AM ( 2 ) 

34470 AM3 ( J ) =AM CIO) 

34430 D0 27 JIL=3,9 

34490 27 AM ( J J L ) = AM ( 2 ) 

34500 IFCJ.EQ.l) G0 T0 28 

34510 G0 T0 29 

34520 31 CONTINUE 

34530 AHZRC J+l )=AHZRC JJ-CC/CCELM 

34540 AZR=AHZRC J+l ) 

34550 PC J+l )»FPRESCAZR,AT) 

34560 AHL0S s DNH*W0RH 

34570 AZR*AHZRC1) 

34530 CEL=CCHR 

34590 AHL*AHZRC J+l ) 

34600 RETURN 

346 1 0 END 

34620C 

34630C 

34640C 

34650 SUBROUTINE BETABNC JJST) 

34660 COMM0N P0, ZAD, EAD, AZR, ALIF, 0GP, GOLD/ TFMX, TFME, CBL, CEL, AHL, 
34670* AKA* AK1 , AK2, WHY, BSRCH, TMNF , BTB0L, BTE0L, SADC, TCCE, SADT, TCT, 

34680* SADH, TCHE, ZADC, ZADT, ZADH, GTH, AGL, CLTH, CGL, Z0P, CCLD, AC0EF , 

34690* AEXP, BUZ# P0W I NT* TIN, DELT 

34700* , TCDGPC I 1 ),TIMC25),TFUELC 1 1 ),GCI 1 > , ATFU ( 25 > , ADCE C 5 > , 

34710* ADT ( 5 ), ADHE C 5 ), TBDC 5 ), AMC 10), TPPKC25 >, PC25 >, BTNAC 1 1 >, 

34720* BTCDGCl 1 ),BDTGPC1 1 >,BTFUEC1 1 ),£0TNC1 1 >,E0TCC1 1 >,£0DTC1 1 ), 
34730* EOTFC 1 1 > t TNAKC 1 1 ),DTGPC 1 1 )> AHZRC 25) 

34740* » H2LC 2 1 it HYLC25 it HRCC C 25 > , TSWELLC 25 > , G 00LC23 ) 

34750* * 0SSW C 25 > , BUSW C 25 ), TINZ C 25 ), DELT ZC 25 ) t P0WZC25 ) 

34760* t AMI (25 it AM2 C 25 ) , AM3 C25 ) 

34770* t TCLADC 25* 1 1 ) 

34730* t KPT.* NJK, JSTP, NSTEPS 

34790 FBETAC Z) a 16l5»97+89» 9747*C ALOG C Z) >+3 . 88 I 0*C AL0G C Z > > **2 

34800* +0*1 3459+C AL0G C Z ) ) **3 

34810 BTBOL=FBETACPC 1 >) 

34820 PZ=PCJSTP) 

34830 BTE 0L*FBETACPZ) 

34840 IFCBTB0L.LE.TPPKC 1 >>G0 T0 5 

34850 6 K* JSTP- 1 

34860 D0 14 1*2, K 

34870 I D= I 

34880 BTP-FBETACPCI)) 

34890 IFCBTP.LE.TPPKC I ) )G0 T0 17 

34900 14 CONTINUE 

34910 26 TFEL-TPPKC JSTP) 

34920 IFCBTE0L.LE.TFEDG0 T0 9 
34930 G0 T0 2 

34940 5 PRINT 22, T IMC 1 it BTB0L, TPPKC 1 ) 
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CONTINUED 


34950 G0 T0 6 

34960 17 PRINT 22,TIMC ID), BTP, TPPKC ID) 

34970 60 T0 26 

34980 9 PRINT 22, TIMC JSTP>,BTE0L, TPPKC JSTP) 

34990 60 T0 2 

35000 22 F0RMATC "CAUTI0N!! BETA PHASE FUEL HAS OCCURRED AT TIME'*, 
3501 0* F8.1," HOURS • " /"BETA FUEL TEMPERATURE IS ",F8.2, /"BUT PEAK 
35020* FUEL TEMPERATURE IS ”,F8.2/> 

35030 2 RETURN 

35040 END 

35050C 

35060C 

35070C 

35080 SUBR0UT I NE XPL 0DE < T I NX, DELTX, TMX, BUM, P0WX, RH0, T I MEND, B0LRH0 
35090* , FPK, UBUKK, SMK, BUM0M1 , BUM0M2, PPBUK 

35 1 00* , FRAC 1 , FRAC2, TDEFK, XEN0NK, HRDK, PDEFK, SMKO 

35 11 0* > 

35120C 

35130 C0MM0N P0, ZAD, EAD, AZR, ALIF, 06P, G 0LD, TFMX, TFME, CBL, CEL, AHL, 

35 1 40* AKA, AK1 , AK2, WHY, BSRCH, TMNF , BTB0L, BTE0L, SADC, TCCE, SADT, TCT, 

35 1 50* SADH, TCHE, ZADC, ZADT, ZADH, GTH, AGL, CLTH, CGL, Z0P, CCLD, AC0EF, 

35160* AEXP, BUZ, P0WINT, TIN, DELT 

35170* , TCDGPC 1 1 ), TIM<25 ), TFUELC 1 1 ),G( 1 1 >, ATFUC25),ADCEC5>, 

35180* ADTC5 ), ADHEC5 >, TBDC 5 ), AM( 10), TPPKC 25 > , PC25 ), BTNAC 1 1 >, 

35190* BTCDGC1 1 >,BDTGPC1 1 >,BTFUEC1 1 ),E0TN(1 1 ),E0TC(1 1 >,E0DTC1 1 ), 
35200* E0TFC1 1 >,TNAK( 1 1 ),DTGP(1 1 >,AHZRC25> 

35210* ,H2L(21 >,HYLC25),HRCCC25>,TSWELLC25), G00LC23) 

35220* , 0SSWC25 ),BUSWC25 ), TINZC25 > , DELTZC25 >, P0WZC25 > 

35230* ,AM1 C25 ),AM2C25), AM3C25) 

35240* , TCLADC25, 11) 

35250* , KPT, NJK, JSTP, NSTEPS 

35260C 

35270 DIMENSION FPC2 1 >, UBUC2 I >, SMC2 1 > , PPBUC2 1 ) , TBARC2 1 >, 

35280* TDEFC21 >,XEN0N(21 >,HRDC21 >,PDEF(21 >, RL0SSC2 1 >, RH0C2 1 ) 
35290* , TMX ( 2 1 > , DELTX (21 ),BUM<21 ),P0WX<21 ),TINX<21 > 

35300* , 0HR ( 2 1 > , XMT C 2 1 > 

353 1 OC 

35320 D0 1 1*1,21 

35330 FP< I )=BUM< I )*FPK 

35340 UBU < I ) *8UM ( I ) *UBUKK 

35350 SMC I )*SMKO*( 1 . -EXPCBUMC I ) /SMK) ) 

35360 PPBUCI > -PPBUK* <F RAC J +EXPC-BUMC I >*21 9.O5*BUM0M1 )* 

35370* FRAC2*EXPC-BUM< I > *2 1 9. O5*BUM0M2 > ) 

35380 TBARC I )*TINX( I )+DELTXC I )/2 

35390 TDEF ( I ) *TBARC I ) *TDEFK 

35400 XEN0NC I )*P0WX< I )*XEN0NK 

35410 HRDC I )*P0WXC I )*HRDK 

35420 PDEF C I )*P0WX( I )*PDEFK 

35430 RL0SSC I )*FPC I )+UBU C I > +SMC I >+H2L C I )+PPBUC I )+TDEF ( I ) + 

35440* XEN0NC I )+HRD( I J+PDEF Cl) 
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REX CONTINUED 


35450 RH0C J I-B0LRH0-RL0SSC J > 

35460 J CONTINUE 
35470 NL0C=O 
35480 VALUE»0. 

35490 LJFE=1 +NSTEPS 
35500C 

35510 00 5 J*1 , 2 1 

35520 J 12=22- J 

35530 OHRC J >=RH0C J12> 

35540 5 XMT Cl) *TMX C 1 1 2 > 

35550 TJMEND*FUNCTI C OHR, XMT, LIFE, VALUE, NL0C, 2 > 

35560 T JMEND^T IMEND/8766 • 

35570C 

35580 JFCKPT.LT. 0) G0 T0 78 

35590 PRINT 77, TIMEND 

35600 77 FORMAT ("E0L C YRS 3X, F 1 0 • 3 ) 

35610 PRINT 20, 

35620 00 2 J * 1 , 2 1 

35630 2 PRINT 10, FPC I ) , UBU ( I ) , SMC I ) , PPBU Cl), H2LC I ) , TDEF Cl), XEN0NC I ), 
35640* HRDC I ), PDEF C J ), RL0SS C J I , RH0 C I > 

35650 10 F0RMATC 1 1F6.2) 

35660 20 FORMAT C 

356704 " FP UBU SM PP0T H2L 

35680* TDEF XENON HRD PDEF RLOSS RH0” ) 

35690 78 CONTINUE 

35700 RETURN 

35710 END 

35720C 

35730C 

35740C 

35750 SUBROUTINE GAGC OPT, JSTEP, GG, GP, BU, DIAFU, AMARG, PHIA) 

35760C CALCULATES GAP AND OTHER FUEL SWELLING CHARACTERISTICS 
35770 COMMON P0, ZAD, EAD, AZR, ALIF, OGP, GOLD, TFMX, TFME, CBL, CEL, AHL, 
35780* AKA, AK1 , AK2, WHY, 3SRCH, TMNF , BTBOL, BTE0L, SADC, TCCE, SADT, TCT, 
35790* SADH, TCHE, ZADC, ZADT, ZADH, GTH, AGL, CLTH, CGL, ZOP, C-CLD, ACOEF, 
35800* AEXP, BUZ, POW I NT, T I N, DELT 

35810* , TCDGPC 1 1 ) , TIMC 25 > , TFUELC 1 1 >, UC 1 I ), ATFUC25 ) , ADCEC 5 ) , 

35820* ADT C 5 ) , ADHE C 5 > , TBDC 5 > , AM C 1 0> , TPPKC 25 ),PC25), BTNAC 1 1 ), 

35830* BTCDGC 1 1 ),BDTGPc 1 1 >,BTFUEC 1 1 ),£.0Twv . i ),c.0 TCC 1 1 ),E0DTC 1 1 ), 

35840* E0TFC1 1 ),TNAKC1 1 >,DTGPC1 1 >,AHZRC25) 

35850* , H2LC 2 1 ) ,HYLC 25 > , HRCC- C 25 > , TSWELLC 25 >, G00L C23 ) 

35860* , OSSW C 25 > , BUSWC25 ),TJNZC25), DELTZC 25 ) , P0WZC 25 ) 

35870* ,AM1 C25), AM2C25), AM3C25) 

35880* , TCLADC25, 1 1 > 

35890* , KPT, NJK, JSTP, NSTEPS 
35900C 

35910 DIMENSION A0FC1 1 >,PHIAC1 1 > 

35920 ALAM S 5000. 

35930 GMX=0. 

35940 0SSWC1)=O 
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REX CONTINUED 


35950 BUSW Cl) =0 

35960 G00LC1 >«GG*1000. 

35970 J-JSTEP 

35980 JK* J+ 1 

35990 00 82 JC=1/ 1 1 

36000 82 GCJC)*GG 

36010 IFC0GP.LT.O>G0 T0 10 

36020 IFCJSTEP.GT.1JG0 T0 20 

36030 D0 22 JA*1/ 1 1 

36040 22 A0FCJA>«O. 

36050 20 IFCTIMC JK> .GT.ALIF JTIMC JKJ*ALIF 
36060 AXPJ=C1 • -EXPC - ( TIMC JKJ-TIMC J ) ) /ALAM ) ) 

36070 00 69 K=l, 1 1 

36080 AT=TFUELCK>+459.7 

36090 A0FF =AC0EF *EXPCAEXP/AT) 

36100 DDV=CAOFF-A0F CK) ) *AXPJ 

36110 IFC0PT.GT.OJG0 T0 7 

36120 IF C K. N£ • 7 >G 0 T0 7 

36130 PRINT/ A0FF/A0F(K) 

36140 7 CONTINUE 

36150 IFCDDV.LE.O. >DDV=0. 

36160 A0F C K) =A0F C KJ+DDV 

36170 ABU=BU*PHIACK>*I .45 

36180 BNUP»3.21 *ABU 

36190 IF CK.EQ * 8 > BUSWC J+l >«BNUP*D1AFU*5 . /2.8 
36200 GR»CBNUP+AOF (K> >*DIAFU*0.0I /2.8 

36210 GCK)*GG-GR/2. 

36220 GRR-GR/2. 

36230 GMX=AMAX1 (GMX/GRR) 

36240 IF C GRR *LT • GG JG0 T0 81 

36250 G C K> = 0 • 

36260 8) CONTINUE 
36270 69 CONTINUE 
36280 T SWELL C J > »TFUELC 8 ) 

36290 0SSWCJK>«A0FC8>*DIAFU*5./2.8 

36300 IFC0PT.GT.OJG0 TO 11 

36310 PRINT 57/ CGCKJ/K*2/ 1 1 > 

36320 57 FORMATC 10F7.4) 

36330 11 CONTINUE 

36340 IF CTIMC JKJ .LT.ALIF )G0 T0 10 

36350 IFCOGP.LT. 1 )G0 T0 10 

36360 RG=GMX+AMARG 
36370 DIF =ABS CGG-RG) 

36380 IFCDIF. LT.O. 00001 >G0 TO 10 

36390 GP-RG 

36400 10 JSTEP*JK 
36410 G0LD=CGG-GMX)*1OOO. 

36420 G00LC JKJ*G0LD 
36430 RETURN 

36440 END 
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REX CONTINUED 


36450 SUBROUTINE FINDGA(G00L#L00Gs TIM#MTT* LUMP* GAPTIM* GAP0) 

3646 OC G00L(LUMP) IS THE X- ARRAY! TIM IS THE Y-ARRAY. SINCE G00L 
36470C IS IN DESCENDING 0RDER, HAVE T0 INVERT IT T0 L00G <ALS0 INVERT 
36480C TIM T0 MIT). GAPTIM IS THE CALCULATED TIME AT WHICH GAP»GAP0. 
36490 REAL G00LCLUMP),L0OGCLUMP),TIM<LUMP>,MIT<LUMP> 

36500 D0 I I « 1 j LUMP 
36510 J=LUMP-I+1 
36520 L00G (J ) *G00L< I ) 

36530 1 MITC J)»TIM< I ) 

36540 NL0C=O 

36550 GAPT IM»FUNCT1 CL00G,MIT,LUMP,GAP0,NL0C*2> 

36560 GAPTIM=GAPTIM/8766. 

36570 RETURN 
36580 END 
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BNDHZR 


40000 SUBROUTINE BNDHZR(BBZR.» AZZ.# ATFUEL* TMMMF ) 

40005 COMMON DUMDUM C 933 ).» IBUM(4> 

40010 FBP(Y)=1 .434E-10*EXPC.01806*Y) 

40020 FPRES(X* AT) *EXP( -8. 8455+88. 8901 *X-78. 8961 *X**2+21 .3731*X**3> 
40030& *EXPC (-12.972+9 • 7707*X-2 . 4984*X**2 ) *1 .0E04/AT) 

40040 ADD=.l 

40050 AZZRaAZZ 

40060 3 CONTINUE 

40070 AT=ATFUEL+459 • 7 

40080 PI =FPR£SCAZZR,AT) 

40090 IF ( AZZR • GT • 1 • 7 )G0 TO 6 

40100 BZR=AZZR+ADD 

401 10 GO TO 7 

40120 6 BZR= AZZR- ADD 

40130 7 P2=F PRES CBZR* AT ) 

40140 A® ( AZZR-BZR ) /( ALOG CPI )-AL0G(P2) ) 

40150 B=AZZR-A*AL0G(P1 ) 

40160 BP=FBP ( TMMMF > 

40170 BHZR=AZZR 

40180 AZZR=A*ALOG(BP)+B 

40190 DIF =BHZR-AZZR 

40200 IFCDIF .LT.O.O)DIF®-DIF 

40210 IFCDIF .LT. 0.001 )G0 TO 4 

40220 ADD=ADD*.2 

40230 GO TO 3 

40240 4 AZZ*AZZR 

40250 RETURN 

40260 END 


INTRP2 


1 1 000 
11010C 
1 1 020C 
1 1030C 
1 1 040C 
1 1 050C 
1 1 060C 
11070C 
1 1 080C 
1 1 090C 
1 1 1 OOC 
1 1 1 1 OC 


FUNCT I ON FUNCT2 ( X 1 , X2, Z, NP 1 , NP2, LP1 , XP1 , XP2, IS.# IC ) 
FEBRUARY 4. 1 972 

TWO-DIMENSIONAL LAG RANG IAN INTERPOLATION ROUTINE 
XI = FIRST ABSCISSA ARRAY 

(MUST BE IN ASCENDING ORDER IN ARRAY) 

X2 = SECOND ABSCISSA ARRAY 

(MUST BE IN ASCENDING ORDER IN ARRAY) 

Z = TWO-DIMENSIONAL ORDINATE ARRAY 
NP1 = LENGTH OF Xl ARRAY 
NP2= LENGTH OF X2 ARRAY 

LP1 = TOTAL FIRST DIMENSIONAL STORAGE FOR Z ARRAY 
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INTRP2 CONTINUED 


11120C XPl = VALUE 0F FIRST ABSCISSA F0R I NTERP0LAT I ON 

11I30C XP2- VALUE 0F SECOND ABSCISSA F0R INTERPOLATION 

1 1 1 40C IS REMEMBERS PLACE IN TABLE (MUST BE SET EQUAL 

I1I50C T0 ZERO BEF0RE FIRST CALL T0 FUNCTI0N) 

11I60C IC C0NTR0LS I NTERP0LATI 0N METHOD* 

I 1 I TOC » 1 F0R LINEAR INTERP0LATI 0N 

11180C a 2 F0R LAG RANG IAN I NTERP0LAT 1 0N 

1IJ90C FUNCT2 SET T0 C 0RRESP0NDING ORDINATE VALUE 

11200C CFUNCT2 WILL EXTRAP0LATE ) 

J1210C 

11220 DIMENSION XI ( 1 ), X2C 1 ), ZCLP1 , NP2 ), YC20) 

11230 DIMENSION EDC 3 ), XA( 3 ) , YAC3 ) 

11240 EQUIVALENCE (EDC 1 ) * XAC 1 ) ) 

11250 NINT-3 

11260 ISJ*0 

11270 IF(NP1.GE.3 .AND. IC.GE.2) G0 T0 15 

11280C LINEAR INTERPOLATION 

11290 NN=1 

11300 1=0 

11310 IF(NPl.LE.l) G0 TO 38 

11320 NINT-2 

11330 15 NN=MAX0(2/NP1 +2-NINT) 

11340 XP=XPI 

11350 NN1 =2 

11360 IFCIS.LE.l .OR. IS.GE.NP1 > GO TO 102 

11370C TEST WHICH DIRECTION TO SEARCH - MUST FIND THE SMALLEST 

11380C X! GREATER THAN XP, BUT WITH I GREATER THAN 1 AND LESS 

J1390C THAN NP1 COR EQUAL TO NP1 FOR LINEAR INTERPOLATION) 

11400 211 NN1 =1 S 

11410 IFCXP.LT. X! (NNl ) ) G0 TO 101 

11420C SEARCH UPWARD IN TABLE 

11430 102 D0 20 DNNDNN 

11440 JFCXP.LE.X1 CD) G0 T0 10 

11450 20 CONTINUE 

11460 I =NN 

11470 GO T0 10 

11480C SEARCH DOWNWARD IN TABLE 

11490 101 DO 21 J=2. NNl 

11500 I =NN1 - J+ 1 

11510 IFCXP.GT.X1 CD) GO TO 11 

11520 21 CONTINUE 

11530 1=1 

11540 11 1=1+1 

11550 10 SPAN=1 .0E-06*XP 

11560 NN«0 

11570 DO 12 J* 1 » NINT 

11580 NNl *1 + J-2 

11590 EDC J)=XP-X1 CNN1 ) 

11600 12 IF CABS (EDC J ) ) • LE. SPAN ) NN-NN1 

11610 IFCNN.EQ.O) GO TO 90 
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JNTRP2 CONTINUED 


11620 
11630 
11640 
11650 
11660 
11670 
11680 
11690 
11700 
11710 
1 1 720 

I I 730 

I I 740 
1 1 750 
11760 
11770 
1 1 780 
12000 
1201 0C 
12020C 
12030C 
12040C 
12050C 
1206 OC 
12070C 
12080C 
12090C 
12100C 
121 IOC 
12120C 
12130C 
12140C 
12150C 
12160C 
12170C 
121S0C 
121 90 
12200 
12210 
12220 
12230C 
12240 
12250 
12260 
12270 
12280 
12290 
12300 
12310 
12320C 


38 00 40 K*1,NP2 
40 YCK)»ZCNN,K> 

FUNCT2=FUNCT 1 CX2, Y, NP2,XP2, JSt* IC) 

G0 T0 80 

90 00 93 J=1,NJNT 
NNl=J+o’-2 
00 92 K=1.,NP2 

92 Y ( K ) = Z ( NN 1 .» K ) 

YAC J)=FUNCT1 (X2jY,NP2jXP2, JS1 > JC) 

93 XA<J)=X1 (NN1 > 

NN 1 *3 

IFCNINT.LT. 3) NN1=2 
IS1»0 

FUNCT2=FUNCTI CXA, YA,NN1,XP, IS1, IC) 

80 JS=J 
RETURN 
END 

FUNCTION FUNCT1 CX, Y, nP, XP J , J S* JC> 

FEBRUARY 4/ 19/* 

THREE POINT LANGRANGJAN INTERP0LATI0N 
0R 

TWO POINT LINEAR INTERPOLATION 
X * ABSCISSA ARRAY 

(MUST BE IN ASCENDING ORDER IN ARRAY) 

Y * 0RDINATE ARRAY 
NP « LENGTH 0F X A Y ARRAYS 
XPJa VALUE OF ABSCISSA FOR INTERPOLATION 
IS REMEMBERS PLACE IN TABLE (MUST BE SET TO 
ZERO BEFORE FIRST CALL TO FUNCTION) 

IC C0NTR0LS INTERPOLATION METHOD* 
a 1 FOR LINEAR INTERPOLATION 
= 2 FOR LAGRANG IAN INTERPOLATION 
FUNCT1 SET T0 CORRESPONDING ORDINATE VALUE 
CFUNCT1 WILL EXTRAPOLATE) 

DIMENSION XCl), YCl), EDO) 

EQUIVALENCE <EDC1),E1), <EDC2),E2), <EDC3),E3) 

NJNT=3 

IF (NP.GE.3 .AND. IC.GE.2) GO TO 15 
LINEAR INTERPOLATION 
NNa 1 
1=0 

IFCNP.LE. 1 ) GO TO 38 
NINT=2 

15 NN=MAX0C2> NP+2-NJNT) 

XP=XPJ 
NN1 *2 

JFCJS.LE.I .OR* IS.GE.NP) GO TO 102 

TEST WHICH DIRECTION TO SEARCH - MUST FIND THE SMALLEST 
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INTRP2 CONTINUED 


12330C X GREATER THAN XP, BUT WITH I GREATER THAN I AND LESS 

123410C THAN NP C0R EQUAL T0 NP F0R LINEAR I NTERP0LAT 1 0N ) 

12350 211 NN1 *1 S 

12360 IFCXP.LT.XCNN1 )) G0 T0 101 

123 7 0C SEARCH UPWARD IN TABLE 

12380 102 D0 20 I*NN1,NN 

12390 IFCXP.LE.XCI>> G0 T0 10 

12400 20 CONTINUE 

12410 I =NN 

12420 G0 T0 10 

12430C SEARCH DOWNWARD IN TABLE 

12440 101 D0 21 J*2.NNI 

12450 I *NN1 - J+ 1 

12460 IF CXP.GT *X( I > > GO T0 11 

12470 21 CONTINUE 

12480 1*1 

12490 11 I*J+1 

12500C INTERPOLATE FOR FUNCT1 IN Y-TABLE CORRESPONDING TO 

1251 OC XP IN X-TABLE USING THREE POINT LAGRANG I AN 

12520C INTERPOLATION SCHEME 

12530 10 SPAN*! .0E-06*XP 

12540 NN=0 

12550 DO 12 «J*1 . NINT 

12560 NN1 *1 +J-2 

12570 EDC J > =XP-XC NN1 ) 

12580 12 IF CABS C EDC J > > «LE • SPAN) NN=NN 1 

12590 IFCNN.EQ.O) GO TO 90 

12600 38 FUNCT1 *YCNN) 

12610 GO TO 80 

12620 90 E 1 2=XC I - 1 ) -XC I > 

12630 IFCNINT.EQ.3) GO TO 360 

12640 FUNCT I *Y CI>-CYCI)-YCI — 1 ))*E2/E12 

12650 GO TO 80 

12660 360 E13*XCI-1 J-XCI+I ) 

12670 E23*X CII-XCI + l) 

12680 36 FUNCT 1 =YC I - 1 ) *E2*E3/CE 1 2*E 13) 

12690& -YCI) *E1*E3/CE12*E23) 

12700* +YCI+1 >*E1 *E2/CE1 3*E23> 

12710 80 IS*I 

12720 RETURN 

12730 END 
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AXDAT 


100 10 * 

110 85; 

120 16 ; 0 ; 

130 .9552; .9352; 

140 ].638;1; 

150 25;2900;98; .933;0; 

160 14750; 1 ; 2 .41 25; 2 . 5; 20800; 

170 .0015; .0015; 0;0;0; 

180 .0015; .0015;0;0;0; 

190 .0015; .0015; 0;0;0; 

200 0 ; 1 00000 ; 0 ; 0 ; 0 ; 

210 0;4383;8766; 13149; 1 7532; 21 91 5; 26298; 30681 ; 

220 35064; 39447; 43830; 482 1 3; 52596; 56979; 6 1362; 65745; 

230 701 28; 745 1 I ; 78894, 83277, 87660, 

240 0,876.6,2191 .5,4383,8766, 13149, 17532,21915,26298, 

250 30681,35064,39447,43830,48213,52596,56979,61362, 

260 65745, 701 28, 745 1 1 , 78894, 83277, 87660, 92043, 96426, 

270 .01, .05,2000, 

280 .00125, .006,2000, 

290 .01, .05,2000, 

300 1120,1 125, 1 134, 1 146, 1161,1177, 1192,1207, 1219, 1228, 1233, 

31 0 1122,1 129, 1 140, 1 154, 1 I 70, 1186, 1202, 1 2 1 5, 1 226, 1 232, 1235, 

320 1 1 38, 1 171 , 1204,1235,1260, I 280, 1 293, 1 296, 1289, 1274, 1251, 

330 1 135, 1 163,1 191, 1218, 1 241 , 1 260, 1 272, 1 277, 1 275, 1 265, 124«, 

340 1122,11 30, 1143,1158,11 74, 1 190, 1206, 121 9, 1228, 1234, 1235, 

350 0, 1 , *004, 

360 0.31,0.659,0.960, 1 .193, 1 .339, 1.389, 1 .339, 1 .193,0.960, 

370 0.659,0.310, 

380 9.047, 

390 10.797, 

400 -.0434, 

410 .02697, 

420 .09352, 

430 3, 

440 .19, 

450 .81, 

460 .00186, 

470 .00163, 

480 .00093, 

490 .00037, 

500 1.04, 
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SI 


100 

35, 

no 

1, 

120 

1, 

130 

.01 , 

140 

1 • , 

150 

1 , 

160 

*014, 

170 

2380, 

180 

.01, 

190 

2380, 

200 

.01, 

210 

1 .0, 

220 

.005, 

230 

.03576. 

240 

.1519, 

250 

.919, 

260 

0, 

270 

• 3 , 

280 

-.04, 

290 

2.4, 

300 

. 0208 , 

310 

1 .0, 

320 

20, 

325 

5, 

330 

.007, 

340 

7.54, 

350 

5.76E0 

360 

-27000 

370 

-1, 

380 

. 1 02, 

390 

.0015, 

400 

0, 

41 0 

1200, 

420 

90, 

430 

1 1 0, 

440 

1, 
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